Thermal conductivity and specific heat of the linear chain cuprate Sr 2 Cu03: 
Evidence for thermal transport via spinons 
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We report measurements of the specific heat and the thermal conductivity of the model Heisenberg 
spin-1/2 chain cuprate Sr2CuC>3 at low temperatures. In addition to a nearly isotropic phonon heat 
transport, we find a quasi one-dimensional excess thermal conductivity along the chain direction, 
most likely associated with spin excitations (spinons). The spinon energy current is limited mainly 
by scattering on defects and phonons. Analyzing the specific heat data, the intrachain magnetic 
exchange J/fcg is estimated to be ~2650 K. 

PACS numbers: 66.70.+f, 75.40.Gb, 74.72.Jt 



There is a considerable theoretical interest in onc- 
dimensional (ID) Heisenberg spin-1/2 systems because 
they exhibit a number of properties that are entirely 
dominated by quantum-mechanical behavior and have 
no analogues in three-dimensional systems. In partic- 
ular, it has been shown that the Heisenberg 5=1/2 
chain represents an integrable system characterized by 
a macroscopic number of conservation laws.u Gvac. im- 
portant conserved quantity is the energy current p a im- 
plying an ideal (infinite) thermal conductivity along the 
chains at nonzero temperatures, if perturbations from 
impurities, phonons, or an interchain coupling, which al- 
ways lead to non-integrable models, are negligible. It is 
an open question, to what extent a real material may 
be regarded as an ideal integrable system. Probably, 
the most obvious evidence for the predicted anomalous 
heat transport is the recent observation of an unusually 
high quasi-lD magnoiijthcrmal conductivity in the series 
(Sr,Ca,La)i4Cu2404i.EKa The structure of these materials 
contains two building blocks with ID character, namely 
Cu02 chains and CU2O3 ladders, both oriented along the 
same direction. Unfortunately, the dimerisation within 
the chains and a non-negligible interchain interaction in 
this system complicate the analysis of the observed fea- 
tures in terms of an integrable model. 

In this work, we have searched for anomalies in the 
thermal transport of S^CuOs, which is often consid- 
ered as the best physical realization of the ID Heisen- 
berg 5=1/2 model. The crystal structure of S^CuOa 
contains-chains formed by CUO4 squares sharing oxygen 
corners.B The chains run along the b axis and, as shown in 
the inset of Fig. |[ the CUO4 squares lie in the ab plane. 
The intrachain exchange interaction between neighbor- 
ing Cu 2+ ions connected via 180° Cu-O-Cu^bxinds, mea- 
sured as J/kr), is between 2150 and 3000 K.ETlI The ratio 
IsbTn/J, where TV is the 3D Neel temperature, is as 
small as 2 x 10 -3 , reflecting an extremely small ratio 
J' I J, J' representing the interchain interaction. 



Our observations indicate an excess thermal conductiv- 
ity along the chain direction, provided by quasi-lD spin 
excitation (spinons) . According to our analysis presented 
below, its magnitude is limited by scattering of spinons 
on defects and phonons. We find no evidence for a mutual 
scattering between spin excitations and hence it seems to 
be absent or at least negligibly small, in agreement with 
theoretical predictions for integrable models. 

The specimens used in these experiments were cut 
from a single crystal that had been grown by the trav- 
eling solvent floating zone method. The details of crys- 
tal growth_and structural characterization are described 
elsewhere. 119 For thermal transport measurements, three 
rectangular-bar-shaped samples of typical dimensions 
2.5 x 1 x 1 mm 3 with the longest dimension parallel to ei- 
ther the a, b or c axis were prepared. Two additional sam- 
ples, #1 and #2, cut from the same piece, were used for 
specific heat measurements. The thermal conductivity 
was measured using a conventional steady-state method 
as described in Ref. |3|. A standard relaxation technique 
was employed for the specific heat measurements. The 
magnetic susceptibility x was measured with a commer- 
cial SQUID-magnetometer. 

Small amounts of excess oxygen are known to be 
present in as-grown crystals of Si^CuOa, giving rise to 
a Curie- Weiss term in thc n tcmperature dependence of 
the magnetic susceptibility^ due to uncompensated Cu 
5=1/2 spins. In order to study the influence of excess 
oxygen, we annealed sample #2 at 870° C for 72 h under 
argon atmosphere, as described in Ref. ||[ From the re- 
sults of our measurements of x(T), the ratio of the num- 
ber of residual spin-1/2 impurities to the total number 
of Cu ions was estimated to be 1.8xl0 -4 for the unan- 
nealed sample #1 and 6xl0 -5 for the annealed sample 
#2. 

The results of the specific heat (C p ) measurements in 
the temperature range between 1.5 and 22 K are pre- 
sented in Fig. |l| as a plot of C p /T versus T 2 . The solid 
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FIG. 1. Specific heat of Sr 2 Cu0 3 as C p /T versus T 2 . The 
solid lines represent the fit to Eq. (|l|). The inset shows the 
data at low temperatures. The value of Tjv denoted by an 
arrow is from Ref. 18. 



lines in Fig. [j] are fits to the data above 4 K using the 
approximation 



C p = jT + f3T 3 + ST 5 



(1) 
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The parameter values are 7=2.12 xlO 
/3=1.359 xl(T 4 J mole" 1 K~ 4 , and 
J mole -1 K~ 6 for sample #1 and 7=2.06 
J mole- 1 K~ 2 , 0=1.531 xl0~ 4 J mole" 1 K~ 4 , and 
5=1.310 xl0~ 7 J mole- 1 K~ 6 for sample #2. The sum 
(3T 3, + ST 5 in Eq. ([!]) is a common low-temperature ap- 
proximation for the lattice specific heat. The fit values 
of the parameter result in values of the Debye temper- 
ature of © D =441±10 K and 424±10 K for samples #1 
and #2, respectively. Since Sr2Cu03 is an insulator, the 
linear in T contribution is not due to itinerant electrons 
but is ascribed to spin degrees of freedom. 

The elementary excitations of a ID Heisenberg spin 
one-half antiferromagnetic system arejpot S=l magnons 
but S—l/2 topological excitations Jl 1 ] now commonly 
called " spinons" . The applicability of the spinon model 
has been demonstrated experimentpJJii-for several 5=1/2 
chain systems, including S^CuOsiiata The .correspond- 
ing specific heat at T -C J/ks is given byEXllj 



C, = 



2Nk\ 
3J 



(2) 



where N is the number of magnetic ions in the system. 
The fit values of the parameter 7 give J/fcs=2620±100 K 
and 2690±100 K for samples #1 and #2, respectively. 
This may be compared with JrikB —2200 K deduced from 
magnetic susceptibility dataQ'EI and a somewhat larger 
value of J//cb=2850 K, as obtained, from the analysis of 
an optical absorption spectrum. LrlHI The cited suscepti- 
bility measurements covered a wide temperature range 
between 5 and 800 K, whereas the absorption spectrum 
presented in Ref. fj] has been recorded at low temperatures 



(32 K). Hence the discrepancy between the values of J 
may be ascribed to its possible decrease with increasing 
temperatureJl3 Our rather large low-temperature values 
of J are compatible with this suggestion. 

The inset of Fig. [I] reveals anomalies of the specific 
heat with onsets below approximately 3.5 K for both 
samples, indicating some sort of phase transition. For 
the annealed sample the anomaly is shifted to lower tem- 
peratures with respect to the peak for the as-grown sam- 
ple. Both anomalies occur at lower temperatures than 
the Neel temperatures T N fovwwf by ^SrH (4.15 <T N < 
6 K) and neutron scatteringc2l (Tjy = 5.4 K) measure- 
ments, respectively. One possible explanation for this 
disagreement is that, besides the transition to an anti- 
ferromagnetically (AFM) ordered state at T/v — 5.4 K, 
not reflected in C P (T), there is another transition at 
T c < 3.5 K. Recently, two subsequent magnetic phase 
transitions at T c \ — 5.0 K and T C 2 — 1.5 K were ob- 
served for SrCu02, containing similar spin-1/2 pchains 
but assembled pairwise in arrays of zigzag chainsB The 
more likely possibility that is consistent with our ob- 
servations is that our anomalies reflect the AFM tran- 
sition reported in Refs. [l9] and but now shifted to 
lower temperatures because of a smaller number of im- 
purities in the samples. It has been shown that nonmag- 
netic impurities interrupting the spin- 1/2 chains enhance 
staggered spin-spin correlations,!^ a common feature 
of various low-dimensional Heisenberg spin systems.E 3 ] 
A convincing manifestation of this feature is the sta- 
bilisation of the long-range AFM order by replacing 
Cu 2+ with non-magnetic ions in the spin-ladder system 
SrCu203 (Ref. |2~i| ) and the spin-Peierls system CuGe03 
(Refs. p5|pq ). Recently, a field-induced staggered mag- 
netisation near impurities was also observed in Sr2Cu03 
by NMR measurements H3 Therefore, it seems quite likely 
that the experimentally observed Neel temperatures are 
always enhanced via the influence of impurities and ex- 
ceed the value that is given by the ratio J'/ J itself. In- 
deed, for Sr2Cu03|-calculations considering only dipolar 
interchain couplinga yield a value of TV as low as 0.028 K. 

The results of the thermal conductivity k along the a, 
b, and c axes are presented in Fig. |[ For both direc- 
tions perpendicular to the chain direction, k(T) shows a 
peak at T max ~ 20 K and a decrease with increasing tem- 
perature, tending to a T^ 1 variation above T > 200 
This behavior is typical for phonon thermal transport S3 
The thermal conductivity along the chain direction, Kb, 
exhibits the same temperature dependence at T < T max 
but obviously not so at higher temperatures. We sug- 
gest that this difference is caused by an additional quasi- 
1D heat transport along the chain direction, provided 
by spin excitations. For insulators the phonon-phonon 
scattering mechanism leads to k oc T~ n (n ~ 1) at 
T > 0£>. For layered structures, such as Sr 2 Cu03, where 
the layers are perpendicular to the a axis, one may expect 
that in this temperature region the ratio between an in- 
plane and the out-of-plane phonon conductivity is larger 
than the anisotropy of k(T) along two different in-plane 
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FIG. 2. Temperature dependence of the thermal conduc- 
tivity of Sr2CuC>3 along the a, b, and c axes. k s is the cal- 
culated spinon thermal conductivity along the b axis. The 
solid line is the estimated sum of spinon and phonon ther- 
mal conductivities assuming that the spinon mean free path 
is equal to the distance between bond-defects (see text). The 
schematic crystal structure is shown in the inset. 



directions H3 Since the difference between k c (in- plane) 
and K a (out-of-plane) is very small, the difference be- 
tween ki, and k c (both in-plane) should even be smaller. 
This argument is, of course, not valid for the temper- 
ature region near and below T max , where the influence 
of sample boundaries and various defects is important 
and, therefore, the behavior of k(T) is, to some extent, 
sample-dependent. Our fitting of k c (T) and n a {T) of 
Sr2CuC>3, employing the Debye model of phonon ther- 
mal conductivity in a similar way as described in Ref. ^, 
has shown that the influence of boundary scattering on 
the thermal conductivity is negligible at T > 50 K. The 
same calculation indicates that an enhanced phonon con- 
tribution to Kb below 50 K may mostly be traced back to 
an enhanced phonon mean free path which is limited by 
boundary and defect scattering. Based on all these argu- 
ments presented above, we make the crucial assumption 
that the phonon thermal conductivity at T > 50 K is 
almost isotropic. 

In order to single-out the heat transport due to spin 
excitations, the phonon thermal conductivity, averaged 
over the a- and c-directions at T > 50 K, was subtracted 
from the experimental data of The resulting spin 
part k s is also plotted in Fig. ^. In a first approxima- 
tion, valid for any system of quasiparticles, the ID ther- 
mal conductivity is given by the simple kinetic expression 
k s = C s v s l s where C s is the specific heat, v s is the veloc- 
ity, and l s is the mean free-path of the spin excitations. 
The velocity of spinons istj v s = Jair/2h, where a is 
the distance between the spins along the chain direction. 
Since T <C J/ks still holds, Eq. (||) for the specific heat 
is valid, and thus the thermal conductivity of spinons is 
given by the simple equation 
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FIG. 3. Spinon mean free path for Sr2Cu03. The solid 
line is a fit to Eq. (EL The solid squares represent the distance 
between two neighboring bond-defects (from Ref. 30) and the 
dashed line is a polynomial fit to these data. 



- N s a^l s T, (3) 
in 

where N s is the number of spins per unit volume. Eq. (Si 
has been shown to also be valid for ID magnon systems.EHI 
We calculated l s using Eq. (||) and taking into account 
small (maximum 6.5% at 300 K) deviations^ of C S (T) 
from linearity. The calculated spinon mean free path is 
shown in Fig. |[ It tends to reach a constant value at low 
temperatures and decreases with increasing temperature. 
Assuming that the different scattering mechanisms act 
independently, the inverse total mean free path of spinons 
l^ 1 may be written as a sum ^Ijl of the inverse mean 
free paths produced by each scattering mechanism. It 
turns out that 

l-^ATexpi-T'ffl+L- 1 , (4) 

with the parameters A=8.2xl0 5 m -1 K _1 , T* = 186 K, 
and L=7.86x 10~ 8 m, reproduces our results in the whole 
temperature region between 50 and 300 K with high ac- 
curacy, as may be seen in Fig. ||. If the first term on 
the right-hand side of Eq. (Q) is to represent Umklapp 
processes, spinon-spinon scattering may be ruled out be- 
cause this case would require T* to be of the order of 
J/k B ~2600K. Since T* is close to 8 D /2=220K, spinon- 
phonon Umklapp processes are the most likely choice for 
this contribution. 

The second and constant term in Eq. (Q) is attributed 
to scattering of spin excitations on defects interrupting 
the Cu-0 chains. In that case the value of the param- 
eter L is a measure for the mean distance between the 
defects. Surprisingly, this distance is much shorter than 
the average distance of 2xl0 4 A between 5=1/2 impu- 
rities, estimated from the Curie- Weiss-type term of the 
magnetic susceptibility. A similar discrepancy has also 
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been encountered in-the interpretation of NMR measure- 
ments on S^CuOsEj Recently, in order to explain some 
peculiar features of NMR spectra of Sr 2 Cu03, Boucher 
and Takigawaa introduced the concept of mobile "bond- 
defects" , which are not related to interstitial excess oxy- 
gen. The calculated mean distance between two neigh- 
boring bond-defects, ..consistent with the NMR data be- 
tween 20 and 60 kJ^3 is shown in Fig. [|. We note a 
reasonably good overlap with our data of l s at T — 60 K. 
The model of Boucher and Takigawa predicts that at 
lower temperatures the interaction between defects is im- 
portant and the number of bond-defects decreases with 
increasing temperature. If our parameter L in Eq. (Q) 
indeed represents the distance between bond-defects, it 
ought to be temperature-dependent below 50 K. Unfortu- 
nately, a quantitative check of this conjecture is difficult 
because of the uncertain subtraction of the phonon back- 
ground at T < 50 K. However, the idea that the bond- 
defects are the main source of spinon scattering at low 
temperatures is, at least qualitatively, consistent with the 
temperature dependence and the anisotropy of the ther- 
mal conductivity also below 50 K. To demonstrate this, 
we calculated the total thermal conductivity (K p h + K s ) 
along the b axis at temperatures between 20 and 60 K 
assuming, first that K p h is isotropic also in this tempera- 
ture range and equal to the average of k c and K a and, sec- 
ond that l s is equal to the distance between neighboring 
bond-defects given in Ref. |32| (the dashed line in Fig. ||). 
The resulting temperature dependence of Kb shown by 
the solid line in-Fig. || is in qualitative agreement with 
the experiment Hj 

Concluding this paper, we return to the question 
whether the present results are relevant vis a vis of intc- 
grable models. We argue that a sizeable quasi-lD ther- 
mal transport mediated by spin excitations does exist 
in Sr2CuC>3. Its magnitude is not exceptional but the 
scatterers, i.e., defects and phonons, limiting the mean 
free path of spin excitations are extrinsic to the magnetic 
system. Our analysis indicates the absence or negligibly 
small influence of spinon-spinon scattering on the ther- 
mal conductivity, in agreement with the predictions made 
for integrable models E Our results imply that a dissipa- 
tionless energy current, expected for systems that fulfill 
the assumptions of the integrable Heisenberg ID 5=1/2 
model, is not robust if perturbations like defects and lat- 
tice excitations interfere. The very high value of the mag- 
netic exchange interaction within the chains and the low 
temperature magnetic phase transition, identified via our 
specific heat measurements, confirm that Sr2Cu03 may, 
nevertheless, be considered as an excellent realization of 
a ID 5=1/2 Heisenberg antiferromagnet. 
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